Spin locking at the apex of nano-scale platinum tips 
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Abstract 

Nanostructures based on platinum, such as small 
clusters or STM-tips, often exhibit an atomistic 
structure that relies upon one or very few strongly 
under-coordinated platinum atoms. Here, we an- 
alyze a paradigmatic example, an apex atom on 
a pyramidal platinum cluster employing the den- 
sity functional theory. We show that such a pris- 
tine platinum tip exhibits a spin polarization of the 
apex atom with a remarkable robustness. Due to 
a depletion of the projected density of states at 
the apex position, the apex-magnetization is effi- 
ciently locked to about 0.6/1b. 

Keywords: Molecular electronics, nano- 
magnetism, platinum electrode, density-functional 
theory 

INTRODUCTION 

Recent years have seen a significant amount of 
cross-fertilization between the research fields of 
spintronics and molecular electronics. And in- 
deed, exploiting spin and charge at the same time 
could open a fascinating venue towards imple- 
menting memory (spin) and logic (charge) on the 
same device. 7 7 For tailoring and controlling spin 
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and charge transport through contacted nano-scale 
devices (molecules, atomic chains) it is necessary, 
if the device properties can be well-defined down 
to the atomic scale. Therefore, a considerable ef- 
fort in Molecular Spintronics is devoted to opti- 
mizing conditions so device properties become (1) 
reproducible and (2) well controlled. 7 The proper 
choice of the contact material is a central part of 
this endeavor, because on the atomic level contacts 
are invasive, in general, and must be considered to 
be a part of the device. 

Here we report a result that is of relevance for 
the use of platinum as a nano-electrode material: 
due to the reduction of the coordination number a 
single platinum atom adsorbed on the bulk mate- 
rial exhibits a very strong tendency towards mag- 
netism. 

As is well known, platinum is an important ma- 
terial from the technological point of view. Be- 
sides being a commonly used electrode in cat- 
alytic chemistry, it has been employed in scanning 
tunneling microscopy measurements. 7 Also, plat- 
inum has been tested as contact material in Molec- 
ular Electronics, especially in break junction ex- 
periments, 7 ? since like gold it forms mono-wires 
when pulling the junction 7 even if bridged by 
small molecules. 7 However quite unlike gold, 
platinum has intriguing magnetic properties that 
may offer qualitatively new possibilities for con- 
tact design. They exist due to the proximity to 
a Stoner instability. While in the bulk form it 
is known to be non-magnetic, platinum develops 
super-paramagnetic behavior in the case of long 
mono- wires 7 and small nanoclusters. 
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Platinum electrodes have been studied by sev- 
eral authors 7 ? ? ? ? in theoretical transport stud- 
ies of nano- wires and break-junctions. Fernandez- 
Rossier et. al ? showed that the finite magne- 
tization of a short Pt chain survives even af- 
ter attaching it to pyramid-like electrodes. It is 
thus clear, that the reduction of the coordina- 
tion number from three-dimensional (12) to one- 
dimensional (2) is sufficient for a magnetic tran- 
sition to take place - at least on the level of 
density-functional-theory, DFT, with conventional 
exchange-correlation functionals. The question, 
that we would like to address here is, whether a 
single apex-atom of such a pyramid by itself would 
be magnetic already. This would be an interme- 
diate case with a coordination number of three 
or four, depending on the specifics of the contact 
pyramid. We are motivated by the observation 
that such geometries are typical for applications 
of platinum as contact material in Molecular Elec- 
tronics. 

Our study suggests that already the single apex- 
atom is magnetic (on DFT-level). Moreover, we 
see that the magnetization invokes a kind of prox- 
imity effect, so that it enters the pyramid down to 
first and second neighbor distances. To the best 
of our knowledge, such a tip magnetization has 
not been reported yet. In the reference 7 the elec- 
trodes were modeled by infinite surfaces, and no 
finite magnetization seems to accumulate in the 
electrode volume. Thiess, et al. 7 used cluster-like 
electrodes in order to simulate a break-junction 
condition, yet the penetration of magnetization to 
the bulk was not reported by the authors. 

Magnetism is treated in our work on DFT-level 
in the all-electron full-potential approach. We 
study platinum pyramids enforcing full- surface 
conditions by applying non-magnetic (i.e. non- 
polarized) boundary conditions on one of the faces 
and varying the total magnetization in the remain- 
ing volume. We observe that the tip magnetization 
decays from the apex towards the non-polarized 
bulk, spanning several atomic sites. The apex 
magnetization is insensitive to the details of the 
boundary conditions, indicating that the magneti- 
zation of the tip stems from the geometric quantum 
confinement. 

Quantum confinement has been discussed as a 
possible agent to drive magnetic instabilities in 



itinerant electron systems before. We mention the 
0.7 anomaly that can show up in the conductance 
of quantum point contacts 7 7 It was explained as 
a consequence of a reduction of screening in the 
contact region due to the confined geometry. The 
loss of screening enhances electron interactions 
and thus can drive magnetic transitions, at least in 
principle. 

Theoretical studies of the phenomenon approx- 
imate the electronic structure by a two dimen- 
sional jellium, 7 which is appropriate for quantum- 
point contacts as realized in semiconductor tech- 
nologies. By contrast, we demonstrate the phe- 
nomenon in a realistic three-dimensional system 
(corner of a Pt crystal) by using atomistic first- 
principles calculations. In the few-atom platinum 
clusters, the origin of spin magnetic moments can 
be attributed to the potential of the core, giving rise 
to localized 3d states with enhanced Coulomb in- 
teraction. Our claim is that the polarization sur- 
vives also in the bulk limit, in regions of geomet- 
ric confinement as tips, corners, necks or surface 
islands and corrugations. 

We emphasize that the basic mechanism, quan- 
tum confinement and breakdown of screening, is 
not restricted to platinum and could apply to other 
metals close to magnetic transitions. So far, an in- 
trinsic electrode magnetization was not reported 
for iridium 7 and palladium junctions, 7 7 mainly 
because electrodes were simulated by plain sur- 
faces. We suspect, that it could exist, never- 
theless, for the physical reasons outlined in this 
manuscript. 

RESULTS 

We have performed extensive density functional 
studies of platinum pyramids of several sizes. 
Here we show results for a pyramid of 18 plat- 
inum atoms, as in [2j All physical quantities are 
calculated in DFT subject to a boundary condition 
of vanishing spin density on the pyramidal face op- 
posite to the apex (see Methods for the details of 
our calculations). 

[1] displays dependency of the total energy on the 
magnetization of the 18-atom Pt-cluster. We find a 
minimum energy for M ~ 4.01 Naturally, this 
corresponds to a magnetization spanning several 
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Figure 1: Total energy of the 18-atom pyramid as 
a function of the total magnetization. Zero energy 
corresponds to the unpolarized solution. 

atomic sites. To inspect the local behavior of spin 
density upon the variation of magnetization, we 
calculate the line density 

Hz) = jj [n t (r)-n ; (r)]dxdj (1) 

where n^^(r) is the spin majority/minority elec- 
tronic density. The z axis goes through the apex 
perpendicular to the interface layer. [2] shows X (z) 
as z goes through the pyramid. Pink areas are re- 
gions of negligible spin-polarization. They are lo- 
cated outside the pyramid, as well as in the inter- 
stitials, giving rise to a ridge-valley picture. This 
reflects the fact that spin polarization is accumu- 
lated in the compact 5 d orbitals, rather than in the 
overlapping sp orbitals. Interestingly, the ridges at 
z — and 2.21 A have constant height independent 
on the total magnetization. The entire variation 
of the magnetization is effectively restricted to the 
third layer (the third ridge starts from red (negative 
A) to end up in green. Finally, the fourth layer is 
the constrained interface. Here only a small arte- 
fact of spin polarization appears. 7 

The ridge-valley structure allows unambiguous 
decomposition of magnetization to atomic layers 
- simply by integrating the line density between 
interstitials. The plot in [3] contains the layer mag- 
netization divided by the number of atoms in the 
layer. Surprisingly, if the spin density is allowed 
to relax freely (i.e. switching off the constraint at 
the interlay er), the highest polarization is located 
in the layer next to the apex. If the constraint is ap- 
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Figure 2: Upper part: Ball-and-stick model* of 
the pyramid. The platinum atoms whose magneti- 
zation is constrained to vanish are darker. Lower 
part: Line spin density X along the z axis of the 18- 
atom pyramid for magnetization varying between 
M=2 and M=7 /ig. Line spin density varies on the 
vertical axis and is color-coded for clarity. Atomic 
layers are positioned at z — 0,2.27, 4.53, 6. 79 A 
starting from apex and visualized by dashed lines 
in the horizontal plane. 

plied, the magnetization decays towards the inter- 
face region as expected, [^demonstrates again that 
the apex is almost insensitive to the interface re- 
gion (i.e. to the constraint). Change in the overall 
constraint is compensated in the deeper layers of 
the pyramid. The robustness of the apex indicates, 
that its magnetization is determined by quantum 
confinement, i.e. geometry of the tip. 

Finally, we demonstrate that the tight locking of 
apex magnetization may be understood as a con- 
sequence of depletion of apex density of states. [4] 
presents the density of states of the 18-atom pyra- 
mid. 7 In the unpolarized case, the region below 
the electro-chemical potential, /I, shows higher 
spectral density, contributed mainly by the 5d or- 
bitals. The chemical potential lies at the spectral 
resonance of a strong apex character. When allow- 
ing for (constrained) spin polarization, this spec- 
tral peak DOS is pushed to lower (higher) energies 
in the majority (minority) spin. At the same time, 
the apex density of states around the chemical po- 
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Figure 3: Magnetic moment per atom in the layers 
0-3 of the pyramid. Zero corresponds to the apex. 

tential is depleted. Hence, small deviations of the 
total magnetization can not be carried over in the 
apex region, since available single-particle states 
have little weight there. 

Discussion 

We discuss the possible impact of two approxima- 
tions underlying our calculations. 

Firstly, we recall that we have to rely upon ap- 
proximate functionals for electronic exchange and 
correlation. Therefore, one expects that our re- 
sults inherit typical artefacts of mean-field treat- 
ments. Such are related to underestimating quan- 
tum fluctuations so that broken-symmetry solu- 
tions are overemphasized. Still, we expect that the 
main features of our calculation survive also in the 
presence of quantum effects. 

Our main argument is that the observed mag- 
netic moment, albeit localized, spans a few atomic 
sites, and has a relatively large magnitude, (M « 
4/1^), which is not easily Kondo- screened by the 
bulk electrons. Hence, one expects that a major 
fraction of the spin will persist down to lowest 
temperatures. 

However, there is a suspicion that this case is 
slightly different from the commonly studied situ- 
ation, where a transition-metal atom sits on top of 
a flat non-magnetic substrate (no pyramid). Mag- 
netic ad-atoms such as Co exhibit Kondo effect 7 
due to anti-ferromagnetic hybridization exchange 
with the paramagnetic host. In the case stud- 
ied here, ferromagnetic-like instabilities at low- 
coordinated platinum atoms require ferromagnetic 
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Figure 4: Total and apex-projected density of 
states of the 18-atom pyramid. The total DOS has 
been divided by the number of atoms. The left 
plot presents an unpolarized solution, the right plot 
presents polarized M = 4.01^. Upper (lower) 
panel shows majority (minority) sector. The ap- 
proximate chemical potential lies at E = —4.89 eV. 
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interactions (like Hund's rule exchange). Hence, 
the specific nature of quantum fluctuations calls 
for further investigations. 7 We remark that the sit- 
uation here bears similarity to the local moment 
formation observed in ferromagnetic atomic con- 
tacts. 7 Iron, cobalt and nickel (1) form localized 
magnetic moments in the contact region, despite 
the itinerant magnetism of these elements in bulk 
(2) and these moments are subject to experimen- 
tally observed Kondoesque fluctuations. 

Second, we have not been accounting for the 
spin-orbit interaction. Indeed, studies of plat- 
inum nano- structures suggest that a significant 
spin-orbit interaction could be present. For in- 
stance, Tosatti and coworkers 7 predicted col- 
losal anisotropy of magnetization of infinite mono- 
wires. 

To address the consequences of spin-orbit inter- 
actions for our results, we recall that the main ef- 
fect of this interaction will be to fix the direction 
of a given magnetic moment. Hence it helps to 
decrease quantum fluctuations and will favor apex 
magnetism working into the direction of our claim. 

CONCLUSIONS 

We have demonstrated that platinum tips develop 
sizable magnetization (4/i#) even in the presence 
of a paramagnetic boundary. The spin density at 
the apex is effectively locked due to quantum con- 
finement, manifesting itself in a local depletion of 
density of states. The magnetization per atom de- 
cays from the apex to the neighboring layers, mim- 
icking a magnetic proximity effect. 

The very presence of magnetization must be kept 
in mind whenever platinum electrodes are used in 
nano- scale transport studies. More generally, our 
results apply to a broader class of nano structures 
where, as we claim, low-coordinated platinum 
atoms drive a local magnetic instability. The unex- 
pected deviation from the normal metal paradigm 
can have consequences important for molecular 
electronics and spintronics. 

METHODS 

We perform first-principles calculations based on 
spin-polarized density-functional theory. 7 The fo- 



cus is on platinum clusters of pyramidal shape 
"cut" from the face-centered crystal of platinum, 
so that the axis points in the crystalline (111) di- 
rection. We use a family of pyramids of 1 1 atoms 
(3 layers), 18 atoms (4 layers, see [2]) and 33 atoms 
(5 layers) to be able to extrapolate finite- size ef- 
fects. The nearest-neighbor spacing is 2.7749 A. 
The threefold symmetry axis (111) of the pyra- 
mid is lifted by including ad-atoms. In order to 
enforce a non-magnetic boundary on the face op- 
posite to the apex (interface layer), the spin densi- 
ties on the atoms of the interface layer are forced 
to vanish. We achieve this by using a variant of the 
constrained density functional theory according to 
the reference. 7 

We add further remarks on technicalities: In 
principle, the constraint specification does not as- 
sure unique ground state spin density and the 
choice of initial density matrix is crucial for the 
convergence of magnetic systems. We have tested 
several initial spin configurations on the atoms and 
did not see evidence for inter-atomic coupling of 
anti-ferromagnetic character. 

To substantiate the impact of basis set size, 
we studied a small 11 -atom pyramid without 
constraints for parameter sets LIGHT/TIER 1, 
LIGHT/TIER2, TIGHT/TIERl, TIGHT/TIER2 of 
the Fhi-Aims package. 7 The total magnetization 
does not differ more than a few percent. Therefore 
we employ the basic basis set light/tier 1 for most 
of the numerical results. Larger basis sets involve 
extended orbitals, these being a drawback in the 
constrained calculations, involving diffusion of 
constrained regions. 

Size effects have minor impact on the accumu- 
lated spin. For the 11, 18, 33-atom pyramids we 
always find M « 4jU# as the energy minimum, al- 
though the distribution of spin over layers changes. 

We confirm the stability of the magnetic solution 
by checking against the GGA+C/ functional 7 (us- 
ing U — 2 eV). In the 11 -atom pyramid, we see 
a small decrease of the apex magnetization from 
0.69 to 0.60 i± B . 
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The spin density on the apex atom of a 
pyramid-shaped platinum tip is insensi- 
tive to changes of the total magnetization 
and stays locked to about 
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